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We have studied the 5=1/2 quasi-one-dimensional antilerromagnet TiOCl using single crystal 
x-ray diffraction and inelastic x-ray scattering techniques. The Ti ions form staggered spin chains 
which dimerize below Tci = 66 K and have an incommensurate lattice distortion between Td and 
Tc2 = 92 K. Based on our measurements of the intensities, wave vectors, and harmonics of the 
incommensurate superlattice peaks, we construct a model for the incommensurate modulation. The 
results are in good agreement with a soliton lattice model, though some quantitative discrepancies 
exist near Tc2- The behavior of the phonons has been studied using inelastic x-ray scattering with 
~ 2 meV energy resolution. For the first time, a zone boundary phonon which softens at the 
spin-Peierls temperature Tsp has been observed. Our results show reasonably good quantitative 
agreement with the Cross-Fisher theory for the phonon dynamics at wave vectors near the zone 
boundary and temperatures near Tsp. However, not all aspects of the data can be described, such 
as the strong overdamping of the soft mode above Tsp. Overall, our results show that TiOCl is 
a good reahzation of a spin-Peierls system, where the phonon softening allows us to identify the 
transition temperature as Tsp = Tc2 ~ 92 K. 



I. INTRODUCTION 

Understanding the properties of low-dimensional spin- 
1/2 systems remains a central challenge in condensed 
matter physics. In quasi-one-dimensional materials, a 
variety of ground-states can be stabilized, depending on 
the subtle interactions beyond the dominant magnetic 
exchange along the spin chain. For example, the low- 
energy physics of the S — 1/2 Heisenberg spin chainfT','^ 
becomes modified once coupling to the three-dimensional 
crystal lattice is taken into account. In the presence of 
spin-phonon coupling, the ground-state is dimerized, cor- 
responding to the spin-Peierls state. [1, 4] If inter-chain 
magnetic interactions are significant, they can promote 
Neel order as an alternate ground state. Even though 
many quasi-one-dimensional 5 = 1/2 spin chain materi- 
als have been studied so far, only a relatively small frac- 
tion are observed to have spin-Peierls ground states. The 
effects of the spin-phonon coupling in these and other low 
dimensional quantum spin systems are important topics 
for further investigation. 

The spin-Peierls transition has been observed in 
many organic charge transfer compounds, such as TTF- 
CuBDT.@, [i, 0, Mm, m In these systems, a static lat- 
tice dimerization has been measured in diffraction exper- 
iments. In addition, x-ray measurements of diffuse scat- 
tering suggested thepresence of a soft phonon associated 
with the transition[6[. However, direct measurements of 
the relevant lattice dynamics are lacking due to the un- 
availability of large crystals for inelastic neutron scat- 
tering. The observation of a spin-Peierls transition in 
the inorganic compound CuGeO^[ll| rekindled interest 



in spin-Pcierls research, as large single crystals could be 
readily grown. However, it is not clear whether CuGeOa 
is a good realization of a conventional spin-Peierls sys- 
tem. The susceptibility does not follow the expected be- 
havior of a spin chain with nearest neighbor coupling pd|. 
and it is believed that significant next nearest neighbor 
coupling along the chain must be taken into account .p^ 
In fact, the ratio of the nearest neighbor to next nearest 
neighbor exchange interaction is close to the critical value 
for spontaneous formation of a spin gap independent of 
a lattice distortion, fl^ Neutron scattering measurements 
also reveal a significant magnetic coupling between neigh- 
boring chains jl3l]. 

Another puzzle in CuGeOs is the lack of an observed 
phonon softening at the transition to the dimerized state. 
Cross and Fisher developed a theory which includes the 
spin-phonon coupling and predicts a softening of a spin- 
Peierls active phonon to zero frequency at Tgp. 4] Neu- 
tron studies were unable to detect the presence of this soft 
phonon, and, in fact, and a hardening of phonon modes 
was observed [3) Recent theoretical analysis sug- 

gests that CuGeOa may fall within the "anti-adiabatic" 
regime (characterized by a large phonon frequency rela- 
tive to the magnetic interaction energ;^!) where a phonon 
hardening would be expected [TtI. fisl. [l9j . The lattice 
dynamics associated with the more conventional "adia- 
batic" regime have yet to be measured in detail. 

Recently, a new 5=1/2 spin chain material has been 
discovered: the transition metal oxide TiOCl. The struc- 
ture of TiOCl is orthorhombic (Pmmn), consisting of TiO 
planes (shown in Fig.[T]) separated by CI layers. The Ti'^^ 
ions are in the electronic configuration with S — 1/2. 
Initially, TiOCl was thought to be an attractive candi- 
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dynamics are studied using high energy-resolution inelas- 
tic x-ray scattering. We observe a soft, zone boundary, 
phonon mode which drives the spin-Peierls transition, 
and we compare our data with the Cross-Fisher theory. 
Finally, Section V contains a discussion and summary of 
our results. 



FIG. 1: Structure of TiOCl at high temperatures in the undis- 
torted state. Left: Schematic of a Ti-O layer showing the Ti 
ion chains along the 6-direction as well as the positions of the 
oxygens. Right: Schematic of the afc-plane which more clearly 
shows the staggered arrangement of adjacent Ti chains. The 
dashed lines highlight the unit cell. 



date for the "resonating valence bond" state !20l|, but it 
has since been found to be primarily a one-dimensional 
magnetic system. Band structure calculations indicate 
that TiOCl is a Mott insulator with the lowest occu- 
pied Ti dxy orbitals pointed toward each other along 
the crystallographic 6-direction. [2l|, |23| The high tem- 
perature magnetic susceptibility is well described by the 
Bonner-Fisher curve, indicatin g a nearest neighbor mag- 
netic exchange of J 



660 K. |21| In the TiOCl crystal 
structure, the spin chains are staggered-every other Ti 
chain is displaced by half a lattice spacing along the b- 
direction. Hence, the inter-chain magnetic interaction is 
frustrated. 

Upon cooling, the susceptibility falls away from the 
Bonner-Fisher curve around T ~ 150 K and exhibits 
anomahes at Tc2 = 92 K and T^i = 66 K.[21] The drop at 
150 K has been attributed to the opening of a pseudo-gap 
[Hi [13 J anomaly at 92 K corresponds to the onset of 
an incommensurately modulated state [H,[2^, 26, 27, 2^, 
and the one at 66 K corresponds to a commensurate 
modulation[29j . Similar behavior has been observed in 
the isostructural compound TiOBr.[33 [M [H [11 At 
first sight, the observed transitions (such as the incom- 
mensurate phase) suggest unconventional spin-Peierls be- 
havior. In order to examine the applicability of a spin- 
Peierls framework, several issues need to be clarified, such 
as the nature of the incommensurate phase and the lat- 
tice dynamics associated with the structural dimeriza- 
tion. 

In this paper, we present x-ray diffraction and inelastic 
x-ray scattering data on single crystal samples of TiOCl. 
As we shall show, the material undergoes a spin-Peierls 
transition in which the structural dimerization occurs in 
the presence of a soft-phonon at the zone boundary. The 
low temperature structural modulations are discussed in 
some detail. The format of the paper is as follows: Sec- 
tion II outlines the crystal growth and thermodynamic 
characterization of the samples. Section III contains a 
discussion of the low temperature structures observed 
with synchrotron x-ray diffraction. Based on measure- 
ments of the superlattice peaks, including higher harmon- 
ics, we develop a model for the incommensurate structure 
described by a soliton lattice. In Section IV, the lattice 



II. SAMPLE GROWTH AND 
CHARACTERIZATION 

The single crystals used in this work were grown by the 
vapor transport method, using the procedure discussed 
previously [21|. The typical size of an individual crys- 
tal was about 50 mm x 50 mm x 20 fim along the a, b 
and c axes respectively. The magnetic susceptibility was 
measured on a sample composed of six single crystals co- 
aligned to within 4 degrees, for a total sample mass of 
about 8 mg. The susceptibility results (shown in Fig.[2Ka) 
and the inset) were performed using a SQUID magne- 
tometer in an applied magnetic field of 2 Tesla. After the 
subtraction of a small isotropic Curie tail, the suscepti- 
bility curves along the three crystallographic directions 
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FIG. 2: Thermodynamic characterization of our sample, (a) 
The magnetic susceptibility upon warming and cooling after 
the subtraction of a small Curie tail. The sharp drops at 
rc2=92 K and Td = 66 K (cooling) indicate the two transition 
temperatures. Inset: Susceptibility measured along the three 
crystallographic directions, (b) The quantity d{xT)/dT (c) 
The specific heat Cp{T). 
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are nearly identical (as shown in the inset), except for a 
small temperature independent offset, most likely due to 
anisotropic crystal field contributions. This behavior fur- 
ther indicates the Hciscnbcrg nature of the spins. More- 
over, the isotropic drop in the susceptibility is a signa- 
ture for spin-singlet formation, as opposed to Neel order. 
Measurements of the susceptibility (Figl2]) upon warming 
and cooling show hysteretic behavior around Td = 66 K, 
which point to a first order phase transition. Hysteresis 
is not observed around the transition at Tc2 = 92 K, in- 
dicating that this transition is second order in nature. 
Below Tci , the susceptibility becomes temperature inde- 
pendent, as seen previously 2l| , consistent with a fully 
formed singlet state. This apparent first order transition 
to the spin-Peierls ground-state seems to differentiate Ti- 
OCl from the other materials in which a second order 
spin-Peierls transition is observed([l H, 0, 11 i, 0, [H ) • 
However, as we will show in section VI, we find that it is 
actually Tc2 = 92 K which defines the spin-Peierls tran- 
sition. 



Also shown in Figl2i;b) is the quantity d{x{T)T)/dT 
which should be proportional to the magnetic contribu- 
tion to the specific heat using the Fisher relation [si', 'sF]. 
Indeed, the measured specific heat Cp(T) in Fig^c) 
shows clear anomalies at Td and Tc2 which find close 
correspondence to anomalies in d{x{T)T)/dT. This sug- 
gests that the origin of the anomalies in Cp{T) have a 
significant magnetic contribution. A rough estimate of 
the magnetic contribution can be obtained by integrat- 
ing Cp/T over T in order to calculate the change in 
entropy over the transition region. At these relatively 
high temperatures, the phonon contribution dominates 
the specific heat and must be subtracted. We fit the 
background in the vicinity of the transition region to 
a fourth order polynomial denoted by the line in the 
Fig. WLc) ■ After background subtraction, the integration 
yields AS* = 0.04(3)/7Vi4rB for the entropy released in 
the vicinity of the two transitions. For comparison, one 
expects S/Nks ~ {2/3){kBT/J) for the entropy of a 
S' = 1/2 uniform chain at low-temperatures (s^]. Using 
J/ks = 660 k[21| and taking T to be an average of 
and Tc2 yields a rough upper bound for the change in 
magnetic entropy of AS/Nks — 0.08. In a previous spe- 
cific heat study, Hemberger et al. [36| estimated a value of 
0.12 for S/Nkb, obtained by integrating over a wider tem- 
perature range. Our estimate is only meant to account 
for the changes in the near vicinity of Td and Tc2, since 
the subtracted background must also contain a weakly 
temperature dependent magnetic contribution. Within 
the errors, our results suggest a significant portion of the 
entropy change at the transitions is magnetic in origin. 




FIG. 3: X-ray diffraction on a single crystal, (a) Longitu- 
dinal scans and (b) transverse scans through the (0, 2.5, 0) 
commensurate position, (c) Scans along H through the in- 
commensurate and commensurate peak positions. The scans 
through the incommensurate peaks were performed at a dif- 
ferent value of K than the commensurate peak. 



III. SINGLE CRYSTAL X-RAY DIFFRACTION 

A. Superlattice Peaks 

Synchrotron x-ray scattering measurements were per- 
formed on our single crystal samples using the X20A 
and X22C beamlines at the National Synchrotron Light 
Source at Brookhaven National Laboratory. The samples 
were cooled using closed-cycle ''^He refrigerators. Upon 
cooling below Td — 66 K, superlattice peaks were ob- 
served at positions displaced from the fundamental Bragg 
peaks by commensurate wavevectors (0, ±i, 0). This in- 
dicates doubling of the unit cell along the 6-direction, 
consistent with a dimerization of the lattice. Scans 
along H and K through the (0, 2.5, 0) superlattice peak 
(shown in Fig. [31(a) and (b)) are resolution-limited, indi- 
cating that the dimerized structure has long-range order. 
The integrated intensities of 45 fundamental superlattice 
peaks were measured and used to determine the dimer- 
ized structure. The intensities were fit using a model with 
two adjustable parameters: S (the Ti ion displacement 
along the 6-direction) and r (the relative shift of neigh- 
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FIG. 4; Temperature dependence of several fitted parame- 
ters for the commensurate and incommensurate superlattice 
peaks, a) The integrated intensity, where the closed symbols 
refer to a commensurate peak and the open symbols refer to 
an incommensurate peak, b) The full-width at half-maximum 
(FWHM) along the ff-direction. Inset: Hysteresis in temper- 
ature of the commensurate peak intensity, c) and d) The com- 
ponents AH and AK of the displacement wave vector from 
the commensurate position. Inset: Normalized displacement 
wave vector for the incommensurate modulation compared 
with a theoretical prediction involving discommensurations. 



boring Ti chains along the 6-direction) . Note that once 
the chains are dimerized, there is no symmetry prevent- 
ing a relative shift r between the two chain sublattices. 
The best fit values are S = 0.03(1)6, and r = 0.05(2)6. 
For simphcity, our model neglected displacements of the 
CI and O atoms, yet we still find good agreement with the 
refinement parameters of Shaz et al. for the dimerized 
state 0. 

The temperature dependence of the integrated inten- 
sity of the (0, 2.5, 0) superlattice peak is shown in 
Fig. m^a) as the closed circles. Upon warming, the in- 
tensity abruptly drops at Td = 66 K, the same tem- 
perature at which our thermodynamic measurements in- 
dicate a first order transition. In fact, the intensity of 
the commensurate superlattice peak also exhibits a ther- 
mal hysteresis, shown in the inset, with a width for the 
hysteresis loop similar to that measured for the suscep- 
tibility. Therefore, both magnetic and structural probes 
indicate the first order nature of transition at Td. 

Upon warming above Td — 66 K, the structure has 



an incommensurate modulation[23|, |25|, |2a, [27|, |23|, as 
shown in Fig. [3{c). That is, between T^i and Tc2 = 
92 K, superlattice peaks are found at positions dis- 
placed from the fundamental positions by wave vectors 
(±Ai?, i ± AiiT, 0). Here, the incommensurate wave 
vectors are shifted from the commensurate wave vector 
(0, i, 0) by the displacement vectors (±A_ff, ±A/C, 0). 
The widths along the _ff-direction of both the commensu- 
rate and incommensurate peaks are plotted in Fig. [H^b) . 
Both sets of peaks are found to be resolution-limited, in- 
dicating long correlation lengths (greater than 2000 A), 
except for narrow temperature ranges in the vicinities 
of Tci and Tc2- Along the if-direction, the peaks re- 
main resolution-limited over the measured temperature 
range. The incommensurate wave vectors continuously 
change as a function of temperature, and upon cooling 
approach the commensurate wave vector. The values for 
ATJ and AK for an incommensurate peak displaced from 
(0, 2.5, 0) are plotted as a function of temperature in 
Figs. Sfc) and (d). The magnitude for AH is about a 
factor of 5 larger than AK at temperature immediately 
below = 92 K. The temperature dependence of the in- 
tegrated intensity for the incommensurate peak is shown 
in Fig. HJa). The solid line denotes a fit to a power law 
(T - Tc2Y with /3 = 0.3. Unlike the transition at T^, 
the transition at Tc2 appears second order in nature. 



B. Model of the Incommensurate Structure 

Since the magnitudes of AH and AK for the incom- 
mensurate phase are much smaller than i (in recipro- 
cal lattice units), the structure may be thought of as an 
incommensurate modulation of the dimerized structure. 
In CuGeOs, the dimerized lattice of the ground state 
of becomes incommensurate when a large enough mag- 
netic field is applied. The incommensurate structure has 
been interpreted as a lattice of solitons which proliferate 
along the chains, separating locally dimerized regions. In 
CuGeOa , the soliton lattice produces a modulation along 
the chain direction only, hence the incommensurate wave 
vector has a non-zero component along the chain direc- 
tion only. In TiOCl, the incommensurate wave vector is 
displaced from the commensurate position by two com- 
ponents: AH and AK, hence, it is slightly more compli- 
cated. 

As an initial comparison, the temperature depen- 
dence of the displacement wave vector magnitude 
AQ = \/Ai/2 + AK'^ is plotted in the inset of 
Fig. O^c) and compared to a Landau theory predic- 
tion. Incommensurate-to-commensurate lock-in transi- 
tions have been observed in charge density wave systems, 
such as TaSe2 37, 38], upon cooling. McMillan ^39,] used a 
Landau expansion of the lattice free energy to calculate 
the temperature dependence of AQ assuming that the 
deviations from the commensurate wave vector take the 
form of long wavelength phase distortions or discommen- 
surations. The line in the inset of Fig. EIc) is the theory 
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prediction [3£ 



T = 5B K T = 66 K 



d{t) = 4.62/[4.61 + /n(l/t)], 



(1) 



where t = {T - Tci)l{Tc2 - Td) and 5{t) = 
AQ{t) / AQ{t — 1). This prediction is a universal func- 
tion of and there are no adjustable parameters in 
the comparison. The theory roughly captures the sin- 
gular behavior as (5 — + near T^i, suggesting that the 
incommensurate phase in TiOCl may be described us- 
ing a model of discommensurations separating commen- 
surate, dimerized regions. Intensity contour plots at 
different temperatures are shown in Fig. El^a), showing 
how the incommensurate peak positions converge to the 
commensurate position as the temperature is cooled to 
Tel. The observed broadening of the peaks just above 
Tci makes it difficult to comment on possible phase 
coexistence [245 or a discontinuous change in wave vector 
at the incommensurate-to-commensurate transition. 

To obtain a more specific picture of the incommensu- 
rate modulation, we measured the intensities of several 
incommensurate superlattice peaks in the vicinity of the 
(0, 1.5, 0) position. At T = 79 K, we find four peaks near 
(0, 1.5, 0) displaced by wave vectors {±AH,±AK,0) 
as depicted in the intensity map of Fig. [5{b). All four 
of these peaks have the same magnitude for the dis- 
placement vector, and we call these the first harmonic 
peaks. We were able to observe third harmonic peaks 
displaced by (±3Ai7, -SAif, 0) from (0, 1.5, 0), also 
shown in Fig. Eljb). This is the first time the higher 
harmonic peaks have been observed, and their positions 
and intensities provide crucial insight into the behavior 
of the incommensurate modulation. No signal above the 
background could be observed at positions displaced by 
{±3 AH, +3AK, 0) indicating that the intensities of those 
peaks are very weak. These observations allow us to rule 
out a model consisting of the superposition of two inde- 
pendent modulations along the a and b directions, as this 
would not yield third harmonic scattering at the observed 
positions. 

We analyze the intensities using a model which de- 
pends on the real-space displacement vectors Ui of the 
atoms in TiOCl from their equilibrium positions. The 
resulting structure factor is 



(2) 



where fi{Q) is the form factor for atom i and the sum- 
mation is performed over an enlarged supercell for the 
modulated structure. For simplicity, we only consider 
displacements of the Ti atoms, leaving the O and CI 
atoms at their equilibrium positions, and relegate the 
displacements to be along the ^-direction. 

To construct the observed in pattern in Fig. [5ljb) , we 
begin by noting that the continuous evolution of AH 
and AK with temperature indicates that the modulation 
is truly incommensurate. Hence, we first consider the 
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FIG. 5: (a) Series of intensity contour plots showing the tem- 
perature evolution of the superlattice peaks (b) Contour plot 
of the incommensurate peak data at T = 79 K showing both 
the first and third harmonics, (c) Calculation of the inten- 
sity based on our model of the incommensurate structure, 
as discussed in the text, (d) Real-space picture of the Ti 
atom positions in the incommensurate structure. Alternate 
Ti chains have different colors (black and red) for easier view- 
ing. The shaded (unshaded) domains delineate approximate 
regions with the same relative shift +t (— r) between adja- 
cent chains. The lightly colored atoms (yellow) on the domain 
walls denote solitons. 



following sinusoidal form for displacements along the b- 
direction for a single chain: 



u[n) ~ cos ( 



(3) 



where the integer n refers to the nth Ti atom along b, 
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5 is the amplitude of the distortion, and A is the wave- 
length of the modulation (in units of b). The case of 
A = 2 corresponds to a simple lattice dimerization and is 
shown in Fig. [U^a) along with a reciprocal space map of 
the corresponding commensurate superlattice peak. An 
incommensurately modulated structure can be obtained 
by making A slightly larger or smaller than 2. Figure [Sl^b) 
shows a modulation with A = 2.04 along with the corre- 
sponding superlattice peak pattern. 

We observe significantly higher intensities for the first 
harmonic peaks displaced by +AK from (0,1.5,0) relative 
to those displaced by —AK (Fig. EKb)). This naturally 
arises when one considers the two different Ti chains in 
the unit cell. By varying the relative phase shift cj) be- 
tween the modulations on the two chains, different in- 
tensity ratios can be obtained. Figure ^c) shows the 
diffraction pattern corresponding to a unit cell with Ti 
chains whose modulations differ by a relative phase shift 
(j) — O.Stt. The resulting structure factor reproduces the 
observed intensity asymmetry. 

Since the modulation is not simply along the b- 
direction (both AH and AK are non-zero in the incom- 
mensurate phase), we add an additional phase difference, 
^, between neighboring pairs of Ti chains. The final ex- 
pression for the displacements is written as 

u{n, m) — 6 cos — ^ i'P + ^fn^ (4) 

where j = 0, 1 labels the two different Ti chains within 
the unit cell and the integer m refers to the mth pair 
of Ti chains along the a-direction. A schematic for this 
modulation is shown in Fig. [IKd). The addition of the 
phase shift has the effect of rotating the modulation 
direction. The calculated third harmonic peak positions 
fall on the same line connecting the first harmonic peaks 
on either side of (0, 1.5, 0), as seen in the data. This is 
direct evidence for a one-dimensional (single wave) mod- 
ulation describing the incommensurate phase of TiOCl. 
Now, there is a degeneracy between displacements with 
opposite signs in front of (f) and ^. This gives rise to two 
twin domains. By considering the other modulated twin 
domain, the pattern in Fig. EJ^c) is obtained. There is 
good agreement with the calculated intensity ratios in 
Fig. Efc) and the observed ones in Fig.[5]Jb). In addition, 
all of the peak positions are reproduced. We can also de- 
duce the population factors of the two twin domains to be 
about 60% and 40%. In order to reproduce the relative 
intensities of the third harmonics, we used a value for 5 
of 0.26, which is unrealistically large. This indicates that 
deviations from a pure sinusoidal modulation exist, as 
we discuss below in the context of the higher harmonics 
peaks. 

Figure [5l^d) shows the real space positions of the Ti 
atoms in our model for the modulation (in a single twin 
domain). The adjacent staggered Ti chains are colored 
differently (red and black) for easier viewing. A local 
dimerization of the lattice is readily seen throughout the 
plotted structure, where adjacent Ti chains appear to 



Distorted Depiction of 

Structure Modulation 




K=2M, (p=0.3ii, i;=0.45 




m=0 1 2 



FIG. 6; Real-space depictions of structural modulation along 
with corresponding reciprocal space maps of the structure 
factor, (a) Sinusoidal modulation for a single chain which 
yields a commensurate structure, (b) Incommensurate mod- 
ulation obtained by changing the wavelength of the sinusoid, 
(c) Modulation for two adjacent Ti chains, where the value 
for the relative phase shift (f) yields an asymmetry in the in- 
tensities of the first harmonic peaks, (d) Model for the incom- 
mensurate modulation in TiOCl, obtained by phase shifting 
the modulation of each unit cell along the a-direction. 

have a small relative shift r along the chain direction. 
For convenience of viewing, each shaded (unshaded) re- 
gion approximately delineates a single phase domain con- 
sisting of the same relative shift -|-r (— r). Note that 
shift T is not a parameter in the above model, but it 
naturally arises once and ^ are specified. On the 
boundaries between domains (discommensurations), we 
find Ti atoms which are nearly undimerized (colored yel- 
low in Fig. El^d)). These correspond to locations where 
u{n, m) ~ and may be referred to as solitons. Due to 
the sinusoidal form of the modulation in our calculation, 
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FIG. 7: (a) and (b) Scans through the commensurate, first 
harmonic incommensurate, and third harmonic incommensu- 
rate peaks along the H and K directions. The center of each 
peak was offset in order to plot the three scans on the same 
axis, (c) Ratio of the harmonic peak intensities /3//1 as a 
function of the displacement wave vector AQ. The line is the 
prediction from the soliton model as described in the text. 



the solitons are not well defined, being spread out over 
many Ti atoms. The other twin domain (with opposite 
signs of and ^ is similar to the structure in Fig. [H^d) 
but with domain walls which tilt in the opposite sense 
with respect to the vertical chain direction. 

Figure [Tl^a) shows scans along H and K through the 
first and third harmonic peaks at T = 79 K, along with 
corresponding scans through the commensurate peak at 
T = 10 K. The center of each peak was offset in order 
to plot the three scans on the same axis. The first har- 
monic peak and the commensurate peak are resolution- 
limited, whereas the third harmonic peak shows a clear 
broadening. In CuGeOa , Christianson et al. [i^ observed 
a slight broadening of the field-induced third harmonic 
peaks and attribute it to the effects of disorder. They 
concluded that the incommensurate phase was not long- 
range ordered in the direction of the broadening, even 
though it had a long enough correlation length to appear 
resolution-limited. We believe the same holds true for 
the incommensurate phase in TiOCl in both the a and h 
directions. 

The higher harmonic intensities can be calculated by 
expanding the modulation in terms of the Fourier com- 
ponents; the ratio of harmonic intensities is given by the 
ratio of the corresponding Fourier coefficients. Here, we 



follow an analysis similar to that used to describe the field 
induced incommensurate modulation in CuGeOs (4ll.[4^. 
In TiOCl, we note that within a single twin domain, 
the modulation may be considered to be one-dimensional 
along the direction given by the displacement wave vec- 
tor (Ai/, AX, 0). When scattering from a periodic one- 
dimensional modulation u(a;), the i*'' harmonic will have 
intensity 41] 



d9 sin {iO) sin [Q u (x)] 



(5) 



where 9 = qx and q is the reduced wave vector along the 
modulation direction. For the case of a sinusoidal func- 
tion, we calculate /3//1 ~ 10~® using a displacement 
amplitude of 0.036. The measured ratio at T = 79 K is 
/3//1 ~ 10^^ and, hence, is not consistent with a pure 
sinusoidal modulation. In the limit of a square- wave func- 
tion, we calculate an expected ratio /3//1 ~ 10~^, which 
is larger but closer to the measured ratio. 

For a better fit, we consider the Jacobi elliptic function 



u{l) = e(-l)'sn ( 



(6) 



where / denotes a lattice site and the elliptic modulus k 
can have values between and 1 (A: = corresponds to 
a sine wave and fc = 1 a square wave). The parameter F 
is the soliton half width, and A/2 is equal to the soliton 
spacing. Then, the ratio /3//1 can be written: 



h 
h 



Y 



r2 + y + 1 



(7) 



where Y = exp {-t^K (Vl - fc^) lK(k)\ , and K is the 
complete elliptic integral of the first kind. The soliton 
spacing, A/2, can be written in terms of fc by 



A 



AQ 



2kK{k)T. 



(8) 



At each temperature, the measured values for /3//1 
and the magnitude of the displacement wave vector AQ 
can be used to determine k and F. For the data in the 
range AQ < 0.125, the calculated values for the soliton 
widths were roughly constant with F — 5(1) A, indicat- 
ing relatively sharp domain walls. By keeping the soliton 
widths fixed at this value, the expected ratio /3//1 can 
be calculated as a function of AQ as shown by the line 
in Fig. [TJc). The calculated line agrees reasonably well 
for AQ < 0.125, lending credence to the soliton lattice 
description for the incommensurate phase of TiOCl. In 
Fig. Ulc), the data points at larger AQ > 0.125 have 
/3 //i ratios which are much smaller than the prediction. 
This means that for these temperatures, which are just 
below the onset temperature Tc2, the soliton width is con- 
siderable larger than F ~ 5 A and the modulation is more 
sinusoidal. For the nearest neighbor Heisenberg model, 
Nakano and Fukuyama predict [43j| F = , where Aq is 
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the magnetic gap. Using J = 660 Kf2l'| and Ao = 430 K 
[23I [4J this yields F ~ 8 A, which represents an ap- 
proximate length scale for the soliton width along the 
chains. The experimental value pertains to the soliton 
width along the incommensurate modulation, which in 
TiOCl is not strictly along the chain direction. 

Our results lead to the following picture for the se- 
quence of structural changes. Just below Tc2, the struc- 
ture has a one-dimensional incommensurate modulation 
of an underlying dimerized state. At first, the envelope 
function for the modulation has a relatively short wave- 
length and is approximately sinusoidal. Upon cooling, 
the wavelength increases and the modulation crosses over 
from being sinusoidal to being better described by a soli- 
ton lattice (with solitons of fixed width). The solitons 
form domain wall boundaries (or discommensurations) 
between locally dimerized domains as shown in Fig.[5jd). 
The domains on either side of the discommensurations 
may be roughly characterized as having opposite signs for 
r, the shift parameter between neighboring chains. The 
soliton spacing (spacing between domain walls) increases 
with decreasing temperature and diverges at Td . At this 
point, the system wants to become uniformly dimerized, 
and a single r must be selected. This selection between 
two degenerate ground states with different symmetries 
may help explain the first order nature of the transition 
at Tr^ ■ The above description applies to each twin 
domain. 

Our model, which focuses only on the Ti displace- 
ments, is consistent with the structural refinements of 
Schonleber et aZ.[2^. However, by measuring the higher 
harmonic content, we have obtained additional details on 
the structure. The above picture also gives insight into 
the reason for the incommensurate modulation. |26l. [28j 
Assuming the lattice energy prefers an equal spacing 
between Ti atoms between chains, then for uniformly 
dimerized chains, this preference cannot be satisfied due 
to the staggered nature of the chains. The sinusoidal 
modulation allows the lattice energy cost to be spread 
out over many unit cells. As the dimerization amplitude 
increases, it becomes energetically favorable to form fully 
dimerized domains at the local level, with the cost in lat- 
tice energy being relegated to domain walls. Eventually 
upon cooling, the magnetic energy gain due to dimeriza- 
tion dominates, and the lattice prefers a uniformly dimer- 
ized state with the value for t which minimizes the lattice 
energy cost. 



IV. INELASTIC X-RAY SCATTERING 

A. Lattice Dynamics 

The lattice dimerization which accompanies the spin- 
Peierls transition is expected to be a displacive struc- 
tural phase transition. Such displacive transitions are 
characterized by a softening of a zone-boundary phonon 
mode which has a polarization similar to that of the 




Energy (meV) 

FIG. 8; Enerj;y transfer scans through the phonon excitations 
at different Q positions along {Q K Q), at temperatures of (a) 
T = 80 K, (b) T = 100 K, and (c) T = 300 K. The solid 
lines are fits to a damped harmonic oscillator structure factor 
convoluted with the instrumental resolution as described in 
the text. 

static distortion. [4^ |43, This mode softens progres- 
sively as the temperature is decreased, until the mode 
energy reaches zero, and the distortion becomes frozen 
in. In their theory for the spin-Peierls transition, Cross 
and Fisher have calculated the expected behavior of such 
a soft mode. 4] However, to date, the lattice dynamics 
of a zone-boundary soft phonon in a spin-Peierls system 
have not been studied. Here, we present a detailed in- 
elastic x-ray scattering characterization of the lattice dy- 
namics in TiOCl. Thus far the only work on the lattice 
dynamics of TiOCl has been Raman and infrared spec- 
troscop y rn easurements of the zone center optical modes 
[45I I49I |5Q | . The inelastic x-ray scattering technique has 
the distinct advantage in being able to probe the phonon 
modes throughout the Brillouin zone. Also, while inelas- 
tic neutron scattering requires relatively large crystals, 
inelastic x-ray scattering can measure very small crystals 
of interesting new correlated electron systems. 

The inelastic x-ray scattering measurements were per- 
formed using the SRF3ID beamline at the Advanced 
Photon Source at Argonne National Laboratory. An 
in-line nested channel-cut silicon monochromator and a 
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FIG. 9: (a) Dispersion relations for the observed phonon 
modes where the undamped frequency ujo, extracted from 
the damped harmonic oscillator fits, is plotted. The lines 
are predictions from a lattice dynamical calculation, (b) The 
damping parameter F for the lowest energy mode. 



bent silicon analyzer in backscattering geometry were 
used to achieve high energy resolution. By scanning the 
energy transfer u; through the elastic scattering from a 
plexiglass sample, the energy resolution can be deter- 
mined. The instrumental resolution function has the fol- 
lowing approximate form 



,2 /In 2 

+ (1-^)- — 
7 V 71- 



1/2 



X exp 



-4 In 2 



(9) 

where 7 is the full- width and half- maximum (FWHM), 
I„ i s the integrated intensity, and r/ is a mixing parameter 
[5l|. We determined the spectrometer energy resolution 
to be 2.3 meV FWHM. The data we present below have 
been corrected for temperature drifts of the monochro- 
mator and are normalized by monitor counts. 

We performed experiments on two samples of TiOCl, 
consisting of crystals coaligned to within 2-3 degrees with 
total thicknesses of ~ 20 /xm and ~ 40 ^m. The incident 
photon energy was 21.3 keV. Representative energy scans 
are shown in Fig. [5] at temperatures of 80 K, 100 K, and 
300 K. For each scan, Q was held constant at various 
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TABLE I: Interaction parameters used in the shell model cal- 
culation, and the resulting equilibrium crystal parameters. 



points along (0, K, 0) between (0, 1.5, 0) and (0, 2, 0). 
The scattered intensity is a direct measure of the dy- 
namic structure factor S{Q, uj) which is related to the dis- 
sipative part of the response function by the fluctuation- 
dissipation theorem. To model the lattice dynamics, we 
used the damped harmonic oscillator response function: 



SiQ,u;) 



2wF 



2p2 



Ha;) + 1] , (10) 



where Uo is the undamped phonon frequency, F is the 
damping constant, and is an amplitude. In the un- 
derdamped case (wo > F) two inelastic peaks are present 
at positive and negative energy transfers. In the over- 
damped case {ojo < F) the scattering has the form of a 
single peak centered about a; = 0. In the extreme over- 
damped limit (lUo <C F) the dynamic structure factor can 
be rewritten as 



SiQ,iu) 



knTl 



TT 7-^ 



(11) 



where 7 = lo'^/T and the high temperature approxima- 
tion of 1 4- 7i(ti') — kBT/hu! has been used. The solid 
lines in Fig. [S] are fits to the damped harmonic oscillator 
cross section convoluted with the instrumental resolution 
function. In some scans, up to three phonon modes can 
be fit. 

Figure [DJa) shows the fitted values for ojo for all of 
the observed phonons as a function of the reduced mo- 
mentum transfer along the (0 K 0) direction. Note 
that in this scattering geometry, the cross section is such 
that only longitudinally polarized phonons are measured. 
The lowest energy mode corresponds to the longitudi- 
nal acoustic branch. The widths of the phonons in this 
branch are plotted in Fig.^b). At T 300 K (weU above 
Tc2 = 92 K), the width of the phonon mode increases as 
Q approaches the zone-boundary. In fact, at the zone- 
boundary, the fit yields T > luq which indicates that the 
phonon is overdamped. The data at T = 100 K (near 
Tc2) show even more dramatic behavior. The width of 
the phonon mode diverges well before the zone-boundary 
is reached, indicating that the mode is strongly over- 
damped. These results are the first indication that the 
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FIG. 10: Calculated cross sections for the longitudinal 
phonons in the undimerized state at various (0 K 0) posi- 
tions, using the eigenvectors and eigenvalues from the lattice 
dynamical calculation. 

longitudinal acoustic phonon at the zone-boundary is a 
soft mode. 

The solid lines in Fig. Ufa) correspond to lattice dy- 
namical calculations for the longitudinal phonon modes 
based on a shell model using the General Utility Lat- 
tice Program [s^]. The calculation uses phenomenolog- 
ical potentials to model the strong short-range interac- 
tions of the neighboring ions [5^. Four Buckingham po- 
tential interactions serve as a repulsive force between 
atoms, and an additional three-body interaction was used 
to constrain an angle of 90° for the Ti— CI— Ti bonding. 
With these potentials, the system settled into an equilib- 
rium structure nearly identical to the experimental unit 
cell [131 (see Table [I]) . A comparison between the cal- 
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FIG. 11: Energy transfer scans at (0, 1.5, 0) for different 
temperatures. The data have been offset for clarity. The 
T = 15 K data were measured at (0, 1.55, 0) to avoid con- 
tamination from the elastic superlattice peak at a; = 0. 



culated and observed optic phonons in Fig. [Slja) shows 
excellent agreement. However, the calculation signifi- 
cantly overestimates the energy for the acoustic branch 
near the zone boundary, another indication of the soft 
mode. We note that our calculation and others 55] indi- 
cate that this mode is degenerate at the zone boundary. 
The energy eigenvalues and eigenvectors resulting from 
the lattice dynamical calculations allow us to calculate 
the expected scattering intensities using the dynamical 
structure factor convoluted with the instrumental reso- 
lution. This calculation is shown in Fig. (TUl The agree- 
ment with the observed intensities that we measure is 
quite good. Again, deviations from the data are appar- 
ent for low energies near the zone-boundary due to the 
presence of the soft mode. 

The temperature dependent behavior of the zone 
boundary phonon is shown in Fig. 1111 The data for 
T = 300 K and T = 200 K were fit to Eq. [TO] (the gen- 
eral cross section for the damped harmonic oscillator). 
However at T = 150 K and below, the fit no longer con- 
verged, and Eg. llll was used for the strongly overdamped 
case. In this case, it is no longer possible to determine 
the values of lOq and F independently. Useful information 
can be extracted by noting that the integrated intensity 
/ is directly proportional to T/w^. Therefore, a plot of 
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T(K) 

FIG. 12: Behavior of the soft phonon at the zone boundary 
position (0, 1.5, 0). (a) The quantity T/Intensity as a function 
of temperature. The solid line follows the form a{T — Tc) 
with Tc = 88 K. (b) The quantity ujI/T as a function of 
temperature. The solid line serves as a guide to the eye. 
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T/I versus T (shown in Fig. fl^ a)) should have the same 
temperature dependence as to^. The solid line has the 
form Lul — a{T — To) which is the result from mean-field 
theory for a soft-phonon transition, [13, 13 where To is 
the structural transition temperature and a is a constant. 
The fitted line yields To — 88(2) K which closely matches 
the temperature Tc2 = 92 K for the onset of the incom- 
mensurate superlattice peaks. The quantity ujo/^ can 
also be extracted and is plotted in Fig. [T2Fb). We see 
that this quantity plunges to zero around 90 K. This is 
consistent with — > (accompanied by an increase in 
r) as expected for a soft phonon transition. Therefore, 
our inelastic x-ray scattering data provide clear evidence 
for a soft phonon transition occurring in TiOCl. Specifi- 
cally, we find that a softening of the longitudinal acoustic 
phonon at the zone-boundary drives the structural tran- 
sition to the incommensurate (nearly dimerized) state at 
Tc2 = 92 K. 



B. Comparison with the Cross-Fisher Theory 



phonon frequency due to the presence of spin-phonon 
coupling or other anharmonic interactions. 

Cross and Fisher (CF) calculated the polarizability 
n(Q,a;) which takes into account the spin-phonon cou- 
pling. The phonon system was treated using a mean-field 
random phase approximation (RPA), whereas the spin 
dynamics were calculated in a non-perturbative fashion, 
more accurate that the Hartree approach. Their treat- 
ment begins by considering a set of non-interacting spin 
chains, where each chain is governed by the nearest neigh- 
bor Heisenberg Hamiltonian: 



E 



J{l,l + l)Si-Si+i. 



(17) 



Here, the planes of atoms perpendicular to the chains are 
constrained to move together. Assuming a linear depen- 
dence of J on lattice distortions, the exchange coupling 
may be written 



Jqlb 



), (18) 



The preceding analysis relied on fits to a phenomeno- 
logical damped harmonic oscillator model. While this 
provides strong evidence for the presence of a soft mode, 
details regarding the microscopic interactions, such as the 
spin-phonon coupling, cannot be obtained. A general re- 
sponse function which includes such interactions has the 
form 



A{ijj) = -Im 
n 



1 



w2-rj2_n(Q,cj) 



(12) 



where fJo is the harmonic or bare phonon frequency, and 
n((5, uj) is the polarizability (or phonon-self energy). The 
dynamic structure factor is then given by 



S{Q,Lj)^\S{Q)\^[l + n{uj)]A{uj). 



(13) 



By comparing this cross section with that for the damped 
harmonic oscillator, one can sec that Il{Q,uj) contains 
information regarding the damping. Since Il{Q, lu) is the 
only term which is complex, it can be broken into real 
and imaginary parts 



A{lo) 



Im 



U{Q,u;) 



(c.2 



Re 



(14) 



n(g,w)J)2 + /m^n(Q,w) 

If we then make the substitutions 

Re [n(Q, u)] - ujHQ, lo) - nl{Q) (15) 
/m [n(Q,w)] =cjr(Q,cj), (16) 

we recover the dam ped harmoni c oscillator response func- 
tion, where ujq = V^o ^ i?e[n]. Here, represents the 
"quasi-harmonic" frequency which differs from the bare 



where q is the reduced wave vector along the chain di- 
rection, g{q) is the spin-phonon coupling, and Q{q) de- 
notes the phonon normal mode coordinates. The spin- 
phonon interaction leads to an expression for the po- 
larizability that depends on the dimer-dimer correlation 
function ([(S; • S;+i)( , (So • Si)j^q]) which CF calculate 
using bosonization. It is also possible to calculate the 
dimer-dimer correlation function using conformal field 
theory |56|, which yields an identical result in the long 
wavelength limit. The CF polarizability is calculated to 
be i 



UcF{q,(^) = 

X/l 

where 



-0.37 1(1 -e^«'')5(g)P 



LU + c{q — 2kf) 
2^ 



h 



c{q~2kf) \ 1_ 
J T' 
(19) 



2ttT 



h{k) 



1 r(i + izfc) 

\iky 



87rri 



(20) 



where 2kf = -n/b and c — 7rJ5/2 [5^ is the spin-wave 
velocity. In this treatment, the condition that defines 
the spin-Peierls phase transition temperature, Tsp, is 



ncF{q = 2kF,Lo^0,T^Tsp)^0. (21) 



This can be solved to get an expression for the spin- 
phonon coupling for the zone boundary mode: 



g — Qo\j —Tsp. 



(22) 



We now compare our data with the results of Cross 
and Fisher by fitting our data with the dynamic struc- 
ture factor that explicitly includes the CF polarizabil- 
ity Hci?(q', w) fEq. I19p convoluted with the instrumental 
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FIG. 13: Energy transfer scans of the low energy phonon 
mode at the zone boundary. The hnes denote fits to the 
dynamical structure factor calculated using the Cross-Fisher 
theory. 



resolution. Only data in the energy range -20 meV< 
Lu < 20 meV was fit, thereby focusing on the longitudinal 
acoustic mode. We first fit data taken the zone bound- 
ary position (0,1.5,0) at several temperatures, as shown 
in Fig. [T31 At the zone boundary, the CF polarizability 
reduces to 



IlcF{q,io) = -0.7 V 



V27rT 



2 1 
T' 



(23) 



Therefore, the only adjustable parameters in the dynam- 
ical structure factor are the spin-phonon coupling g and 
the bare phonon frequency fio- Since scales with the 
CF polarizability, a non-zero g will affect both the peak 
width and the shift of loq from Qq. Away from the zone 
boundary, the polarizability also depends on J, therefore 
J and g cannot be fit independently in general. 

The solid lines in Fig. [13] represent fits to the dynam- 
ical structure factor at q = 0.5 based on the CF theory. 
Here, we let g vary with temperature, but kept ilp fixed 
at the best fit value of flo = 27(1) meV. Letting g float 
with temperature was necessary in order to get reason- 
able fits, and we will comment on this later. Once g 
was determined for a given temperature, the scans at 
other q positions were then fit with JIq and J as the ad- 
justable parameters. For the data sets at T = 100 K and 
T = 300 K, the analysis yielded a roughly consistent mag- 
netic exchange of J ~ 200 K. This value for J is about 
three times smaller than that deduced from the mag- 
netic susceptibility results j21i]. However, our value for J 
is remarkably close considering that it resulted from mea- 
surements of only the phonon positions and lineshapes. 

The resulting values for the bare phonon frequency Qq 
are plotted in Fig.fWa) as a function of the reduced wave 
vector (0, q, 0). The solid line denotes the prediction from 
the lattice dynamical calculation discussed in the previ- 
ous section. Clearly, the bare phonon frequency ^o{q) 
follows the calculated dispersion curve much better than 
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T{K) T(K) 

FIG. 14: (a) The bare phonon frequency Qq of the longitudi- 
nal acoustic mode extracted from the fits to the Cross-Fisher 
theory. The solid line denotes the shell model calculation, (b) 
The fitted spin-phonon coupling g, normalized to the value of 
g calculated using Eq. 1221 c) The quasi-harmonic frequency 
ujo of the zone boundary mode, extracted from fits to the 
Cross-Fisher theory. The solid line denotes the power-law 
form AVT - Tc with Tc = 98 K. 



the mode energies lJo from the damped harmonic oscilla- 
tor fits (Fig.[T2l^a)). Thus, it appears that the CF theory 
captures the anharmonic effects of the spin-phonon cou- 
pling in TiOCl reasonably well. However, in order to ob- 
tain good fits at all temperatures, g was allowed to vary. 
These values are plotted in Fig. [TW b) , normalized by the 
value of g calculated using Eq. [21] (with VIq — 27 meV and 
Tc2 = 92 K), which yields gcross-Fisher = 75 meV^/^. 
The figure shows that the fitted values of g approach 
gcross-Fisher as T — + Tc2 . Sincc the fits to the dynam- 
ical structure factor were performed without knowledge 
of rc2, this result supports the self-consistency and valid- 
ity of the CF expressions describing the soft phonon near 
Tsp. We also have enough information to determine the 
quasi-harmonic frequency uJo at q = 0.5, which is plot- 
ted in Fig. [TiFc) . We see that the spin-phonon coupling 
reduces lOq from the bare value fio by about a factor of 
2. Also, at the spin-Peierls transition, ujo should fall to 
zero. Indeed, the solid line through the data is a fit to 
the form A^/T — Tc where Tc = 98 K, which is close to 
Tc2 = 92 K. Hence, we identify Tc2 as the spin-Peierls 
transition temperature in TiOCl. The congruence with 
the basic features of the CF theory allows us to conclude 
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that that TiOCl is, in fact, a good reahzation of a spin- 
Peierls system. 

While the CF theory can successfully describe many 
aspects of the phonon data, there are some discrepan- 
cies. Such discrepancies are not completely unexpected, 
since the theory is based on a mean-field RPA treatment 
of the phonons in which fluctuation corrections due to 
the phonon dynamics are absent. As noted above, the 
best fit g is temperature dependent and the fitted value 
of J is a factor of 3 lower than that estimated from the 
magnetic susceptibility. This likely stems from limits to 
the applicable range for HcfC^jw), which is most accu- 
rate in the vicinity of q « 0.5 and for temperatures near 
Tsp- In order to more fully describe the data, additional 
anharmonic effects which broaden the phonon lineshapes 
away from q = 0.5 must be included. In addition, we ob- 
serve the zone boundary phonon remains strongly over- 
damped even at temperatures as high as T = 300 K 
Tgp). These unusual features of the phonon dynam- 
ics of TiOCl require better theoretical understanding. 

V. DISCUSSION AND SUMMARY 

The TiOCl compound has many ingredients which 
make it a particularly ideal spin-Peierls system. First, 
it is composed of weakly interacting 5 = 1/2 spin chains 
which are well described by a dominant nearest-neighbor 
interaction along the chain, [^l*! The weak magnetic cou- 
pling between chains is facilitated by the staggered ar- 
rangement of adjacent Ti chains. Indeed, a lattice dimer- 
ization is observed at low temperatures, concomitant 
with spin-gap behavior observed in the susceptibility. 
However, as a consequence of the staggered Ti chains, an 
incommensurately modulated phase develops prior to the 
commensurate dimerization upon cooling. Our detailed 
x-ray scattering studies have led to considerable new in- 
sight into this incommensurate phase. We have devel- 
oped a model which describes how the Ti chains are mod- 
ulated relative to each other and shows good agreement 
with the data. Our analysis of the higher harmonic scat- 
tering indicates a cross-over from sinusoidal to a solition- 
like modulation as the temperature is reduced below 
Tc2 = 92 K. The incommensurate-to-commensurate tran- 
sition at Tci = 66 K is characterized by a divergence of 
the domain wall spacing and selection of a single dimer- 
ized domain. 

TiOCl is also of interest since it is an inorganic mate- 
rial for which single crystals samples can be grown. Using 
single crystal samples, we have performed high-resolution 
inelastic x-ray scattering measurements of the lattice dy- 
namics. This is a noteworthy example in which x-rays 
rather than neutrons can provide the best information 
on the detailed lattice dynamics of a correlated electron 
system. We have discovered a soft phonon mode which 
drives the spin-Peierls transition. This stands in marked 
contrast to the much studied CuGeOs system, for which 



a soft phonon has not been observed. In TiOCl, the lat- 
tice dynamics suggest that the spin-Peierls temperature 
is associated with the transition temperature Tc2 — 92 K, 
and not Td = 66 K as has been previously speculated. 

Our measurements on TiOCl allow to make an un- 
precedently detailed comparison with the Cross-Fisher 
theory. We find that the calculated polarizability can 
successfully describe the measured phonon cross sections. 
That is, the extracted values for the bare phonon frequen- 
cies and spin-phonon coupling give a consistent picture 
of both the phonon dispersion and the phonon softening 
at Tsp. However, several discrepancies exists, such as 
the small value of J and a temperature dependent spin- 
phonon coupling constant. This points to the necessity of 
adding additional anharmonic effects to fully describe the 
lattice dynamics. Interestingly, the soft phonon mode re- 
mains strongly overdamped at temperatures much higher 
than the observed spin-Peierls transition temperature. 

The presence of the soft mode transition indicates that 
TiOCl falls within the adiabatic regime (i.e., small fio). 
Indeed, the fitted bare phonon frequency of the Peierls- 
active mode, ilo — 27 meV, is smaller than the magnetic 
energy scale J ~ 57 meV in the system. Recent theories 
have expanded upon the Cross-Fisher result and shed fur- 
ther light on the conditions separa ting the adiabatic and 
anti-adiabatic regimes, [l^ [13, Il9l | Gros and Werner 
have shown that within RPA a soft phonon occurs only if 
^o/Tsp < 2.2.0 For TiOCl, we calculate no/Tc ~ 2.9 
(using Tc2 = 92 K as the critical temperature). Hence, 
their prediction suggests that TiOCl falls outside of the 
adiabatic regime. However, a more recent theory by Do- 
bry et al. [58*] goes beyond RPA and takes into account the 
dynamics of the transverse phonons arising from inter- 
chain interactions. They find that soft phonon behavior 
may occur for values as high as ^Iq/Tc ~ 3, provided that 
ijJ±/Tsp for the transverse phonon is not very large [ssj. 
Our shell model calculations indicate that the latter con- 
dition is satisfied in TiOCl. Of course, further experi- 
mental studies of the transverse phonon dynamics would 
certainly be useful. Also, more theoretical work based 
explicitly on the structure of staggered Ti chains would 
be necessary for additional detailed comparisons. With 
its wealth of observed phenomena and relatively simple 
crystal structure, TiOCl is an ideal system for quantita- 
tive tests of theories for coupled spin and lattice degrees 
of freedom in one-dimensional magnets. 
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